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Abstract 



(N 

^^ I We investigate one-matrix correlation functions for finite SU{N) Yang- Mills 

Qs^ I integrals with and without super symmetry. We propose novel convergence 

conditions for these correlators which we determine from the one-loop per- 
turbative effective action. These conditions are found to agree with non- 
Mh. perturbative Monte Carlo calculations for various gauge groups and dimen- 

r~| ! sions. Our results yield important insights into the eigenvalue distributions 

p(A) of these random matrix models. For the bosonic models, we find that 
IT^ ' the spectral densities p{X) possess moments of all orders as A^ — > 00. In 

j_j ' the supersymmetric case, p{X) is a wide distribution with an A— independent 

^ asymptotic behavior p{X) r^ A^^,A^^,A^^^ for dimensions D = 4,6,10, re- 

spectively. 
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Recently there has been renewed interest in dimensional reductions of SU{N) Yang-Mills 
theories. It has been argued that these "toy models", apart from some subtleties, allow to 
recover the full unreduced theories in the large N limit [^ . The supersymmetric reductions 
are relevant to D-brane physics and have also been used in various attempts to define non- 
perturbative formulations of quantum gravity, supermembranes and superstrings. For the 
complete reduction, all space-time dependence is eliminated from the gauge "fields", and 
the continuum Yang-Mills path integral becomes an ordinary multi-matrix integral. The 
integrals appear to be ill-defined due to the flat directions (i.e. commuting matrices) in the 
action. Recent calculations have however uncovered that in many cases of interest these 
matrix integrals converge and thus do not need to be regulated. This was analytically found 
in the supersymmetric case for SU{2) in 0|. In |^, 0| we found by Monte Carlo methods 
that the absolute convergence persists for larger values of N . Furthermore, we numerically 
established (and proved analytically for SU{2)) the convergence properties of the bosonic 
(non-supersymmetric) integrals. Perturbative one-loop estimates in favor of convergence 
were presented for the supersymmetric case in [^ and, recently, in 0, for the bosonic case. 
Mathematically rigorous proofs for A^ > 2 are however still missing. For some applications 
of these integrals see e.g. 0, 0, ||, [|. 

The integrals in question are 
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They correspond to fully reduced D-dimensional Euclidean SU{N) Yang-Mills theory. Here 
A/" is the number of real supersymmetries, and for A/" > the only possible dimensions are 
D = 3,4, 6, 10 corresponding to A/" = 2, 4, 8, 16, respectively. If A/" = (no supersymmetry) , 
a priori all D > 2 are possible and we simply omit the terms with Grassmann fields both 
from the measure and action of eq.(P. For the detailed notation we refer to [§], @]. 

In view of [0], [Q, the necessary and sufficient conditions of existence for the integrals 
eq.(|I|) then appear to be 



D = 4,6,10 and X > 2 } for A/" > 
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As already mentioned in 0], it is only the bosonic D = 3 SU{3) integral which escapes a clear- 
cut classification by numerical means as we cannot exclude almost marginal convergence. 
We have now checked that the susy D = 3 integral is not absolutely convergent^ for any N, 
disproving a conjecture made in |Q. 



^ This case, with AA = 2, is special in that, at least for even A^, the integral is formally zero 



For the large A^ limit, eq.(0) implies existence of the integrals for D > 3 in the bosonic 
case and for D = 4,6, 10 in the supersymmetric case. 

Analytic motivation of eq.(^ for general N can be obtained from the one- loop perturba- 
tive calculations of the effective action as presented in P] , : Decompose the matrices into 
diagonal D{X^) (with matrix elements D{X^)ij = SijX^) and off-diagonal components y^ 
(i.e. (Yfj^ii = 0): X^ = D{X^) +y^. Fix a background gauge, expand the action to quadratic 
order in Y"^ and work out the effective action for the diagonal elements x^ . The result for 
the supersymmetric case reads 
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Here the sum is over all possible maximal (i.e. there are A^— 1 links in G) trees G connecting 
the N coordinates. As indicated in eq.(^), there are actually further terms present to one- 
loop order for D = 6, 10, but they are irrelevant for the powercounting arguments below 
(see [El). For the bosonic case M one finds 
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The one-loop approximation is reasonable for well-separated diagonal components x* (that 
is the "infrared" regime). It is now easily verified that superficial powercounting of all 
integrals on hyperplanes of the integration space in eqs.(^,(§) results, at large x^, precisely 
in the convergence conditions eq.(^. 

Clearly, this perturbative argument is not a proof of eq.(§): First, the cumulative effect 
of all corrections to this one-loop effective action is beyond control. Second, there are 
potentially dangerous regions in the integration space where some of the differences {x^ — x^Y 
are small (thus rendering the one-loop approximation invalid) while others are large (leading 
to a potential divergence). 

An analytic calculation of the D = 4, 6, 10 susy partition functions 2j^j^ has been re- 
ported in [|I^. In this work, light-cone variables = Xi + iX^, (p = Xi — iX^ are used. 
Subsequently, and are treated as independent hermitian matrices (i.e. one actually has 
(p = Xi — Xq, (j) = Xi + Xq), and the integrals computed in |]TD|] correspond to eq.(|I]) 



with Minkowski (as opposed to Euclidean) metric. However, the Minkowski integrals are 
divergent and need to be regulated. This calculation might be rendered rigorous if one suc- 
ceeded in deriving the imposed pole prescriptions directly from the matrix Wick rotation 
Xd ^ tXo. 

For applications, it is important to understand the statistical eigenvalue distribution of 
these random matrix models. The simplest quantity is the distribution for the eigenvalues 
Xi of just one matrix, say, Xi, in the background of the other matrices X2, . . . , X^ 



y]. The integral is however not absolutely convergent, as we have checked numerically with the 
methods used in this paper. This divergence is also present in the analytically tractable SU{2) 
case. Adding a Chern-Simons term, as proposed in |1C], does not improve convergence and the 
D = 3 susy integral modified by such a term is divergent for all N. 
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Here, the average <> is with respect to eq.(|T]). The easiest way to investigate the density 
is to consider the one-matrix correlator^ 
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In an ordinary hermitian matrix model, say with weight exp(— TrXf), the density eq.(^ 
falls off at infinity as p(A) ~ exp(— A^). Therefore all moments eq.(|^) exist. 

In the Yang-Mills ensembles for SU{2), we obtained analytically (along the lines of [@) 
that the correlators eq.(^) exist only for low values of k; we find k < D — 3 (A/" > 0) and 
k < ^ — 2 (J\f = 0). For general N, we obtain the convergence conditions in the one- 
loop approximations eqs.(^,(^. Some care has to be exercised, since the most divergent 
contribution is not obtained by simply counting the overall dimensionality of the integral. 
In fact, a more dangerous infrared configuration stems from only a single (say, the i-th) 
eigenvalue x* becoming large in the D-dimensional space. The only exception being, once 
again, the bosonic D = 3 SU{3) integral. Therefore, apart from this one case, we find 

^ r k<D-3 for ^f>0 

-TrXf)<oo iff (7) 

^^ [ A; < N{D - 2) - f £) + 2 for A/" = 

This means that no moments exist for the D = 4 susy integral, and only the first two and 
six even moments for, respectively, the D = 6 and D = 10 susy integrals! 

We have been able to directly verify our conjecture eq.(0) with the Monte Carlo approach 
of 0], [Q. It may seem straightforward to obtain a direct numerical estimate of the eigenvalue 
spectrum of Xi by generating ensembles of X = (Xi, X2, . . . , Xd) with the statistical weight 



■.^f 



^dA^) = 6^P 



-Tr [X^, X^] [X^, Xy 



Vd,n{X) (8) 



where Vd,n{X) is the Pfaffian polynomial coming from the integration over the fermionic 
degrees of freedom (cf ||^, for the bosonic integrals, this term is simply dropped). The 
eigenvalue spectrum of X can in principle be sampled and its histogram generated. In 
practice, it is however impossible to make firm predictions on the tails of the histogram, as 
they comprise an exceedingly small number of samples. For this reason, we rather generate 
Markov chains of configurations X{t) with the statistical weight 

<'^W = lETrXf4^(X)| (9) 



^ One might attempt to compute such correlators within the powerful framework of ||T^. The 
natural candidate would be (jj^Tvcp^) (where cj) = Xi + iXn). Unfortunately, these correlators 
vanish for all k in the Euclidean integral eq.(||) due to the S0{2) symmetry in the {Xi,Xd) plane 
(M. Douglas, private communication). 



This means in particular that we perform one independent run for every value of k. The 
critical information on the convergence condition is contained in the autocorrelation function 
of our Markov-chain vector X{t) 
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for a fixed value of A (as in [^, we have compactified all the integrands). An absolutely 
non-integrable singularity of the integrand corresponds to an infinite statistical weight con- 
centrated in a region of negligible extension. This simply means that the simulation should 
get stuck, and the autocorrelation function approach 1. An integrable singularity is not 
picked up by the present method. 

We have computed the autocorrelation function for a large number of cases. Among the 
supersymmetric cases, we considered A^ < 6 for D = 4, A^ < 5 for D = 6 and A^ = 2, 3, 4 for 
D = 10. All the analytically known results for SU{2) were easily reproduced. Convergence 
of bosonic integrals was checked for A^ < 9 and D = 3,4, 6. 

As an example, we show in Fig. 1 the autocorrelation functions of typical runs for the 
bosonic integrals with D = 3, N = 6 for k = 2, k = 3 and k = 4. The criterion eq.(|^) leads 
us to expect convergence for k = 2 and A; = 3, as the smallest diverging power is kcru = 3.5. 
This behavior is clearly reproduced in Fig. 1. Notice that the "almost" divergent case k = 3 
gets stuck for very long periods of Monte-Carlo time, whereas the the almost "convergent" 
integral k = 4 settles to a unit value of // only after a very long transient. 
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Fig. 1 Autocorrelation functions ff" versus Monte Carlo time t (A = 1000^) for the bosonic 
integral with N = 6, D = 3, and, from the left, k = 2,3, 4. 



In Fig. 2, we show corresponding plots for the supersymmetric case with D = 6, N = 4. 
Here, we expect the k = 2 integral to be convergent, as clearly found, while k = 4 is 
divergent, as expected. In the marginally divergent case, we have observed the typical 
alternations between "stuck" and "mixing" behavior. 
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Fig. 2 Autocorrelation functions f^ versus Monte Carlo time t (A = lOOj for the super- 
symmetric integral with D = 6, N = A, and, from the left, k = 2,3, 4. 



Notice that the middle plots of Figs 1 and 2 do not appear fundamentally different, even 
though we expect the integrals to behave differently. Our Monte Carlo procedure leaves 
a margin of error which we estimate to be Ak ~ ±1/2. A similar ambiguity affects the 
abovementioned bosonic integral for D = 3 and A^ = 3. 

The conditions eqs.(|^) mean that in the Yang-Mills case the asymptotic behavior of the 
one-matrix eigenvalue distribution eq.(|3) decays algebraically: 



p(A) ~ |A| 



for 



A -^ ±oo. 
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Assuming this to be exactly true we can extract the power a in eq.(|TID: 
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Most strikingly, the decay of the densities in the susy cases D = 4, 6, 10 (p(A) ~ 
A~^, A~'^, A~^^) is independent of A^. It means that the eigenvalue distribution are wide 
even in the N ^ oo limit! This is a most unusual effect for a random matrix model. Ev- 
idently, supersymmetry is responsible for this behavior, as can e.g. be seen by comparing 
the one-loop effective actions for the diagonal matrix elements in the susy (c/eq.@) and 
non-susy (c/eq.(§)) case: In the latter, this effective action is (9(A^^) while in the former it 
is 0{N). 

The bosonic case is much more conventional in that the density becomes concentrated 
in a finite interval at large A^. In this context it is interesting to mention the result of a 
numerical study of the bosonic Yang-Mills integrals which is complementary to the present 
work: In 0, Hotta et.al. investigate the scaling behavior of (-^TrXf) (i.e. eq.(||) with 

k = 1) with A^. They find that, if the variables in eq.([^) are rescaled as X^ -^ N^X^ 
such that an explicit factor of N appears in front of the action, this correlator tends to a 
constant as N —>■ oo. In consequence, the usual 't Hooft scaling of large N matrix models 
is obeyed, and the edge of the eigenvalue support tends to a constant. Nevertheless, the 
observed asymptotic behavior (p(A) ~ A^^^'''^"^-' as N ^ oo) is different from the universal 
exponential decay law of Wigner-type random systems. 
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To summarize, we have demonstrated how statistical information on the eigenvalue dis- 
tributions of Yang-Mills integrals (which are the simplest examples for the so-called "new" 
matrix models) can be obtained. We observed an interesting difference in the tails of the 
spectral distribution between the susy and non-susy ensembles. While suggestive - e.g. we 
find it interesting that the susy eigenvalue distribution tends to stretch out much farther 
- the present study clearly does not yet address such important issues as the full nonper- 
turbative effective action for the diagonal elements of all D matrices or the problem of 
level spacing statistics. These questions will need to be addressed if one intends to apply 
Yang-Mills integrals to string theory and large N gauge theory. 
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